Non dissipative decoherence of Rabi oscillations 
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We present a simple theoretical description of two recent 
experiments where damping of Rabi oscillations, which cannot 
be attributed to dissipative decoherence, has been observed. 
This is obtained considering the evolution time or the Hamil- 
tonian as random variables and then averaging the usual 
unitary evolution on a properly derived, model-independent, 
probability distribution. 



Even though decoherence is a very general phe- 
nomenon p|, it is very difficult to verify it experimen- 
tally because most often the physical nature of the en- 
vironmental degrees of freedom responsible for the deco- 
herence process remains unknown. The only controlled 
experimental verification of decoherence has been given 
by the experiment of Ref. H], in which the progressive 
transformation of a linear superposition of two coherent 
states of a microwave cavity mode into the correspond- 
ing statistical mixture has been monitored. In this case, 
the environmental decoherence has been checked with no 
fitting parameters because its physical origin, i.e., pho- 
ton leakage out of the cavity, was easily recognizable and 
measurable. In this case, it is even possible to control de- 
coherence, i.e., to considerably suppress its effects, for ex- 
ample by using appropriately designed feedback schemes 

i- 

In some cases however, the mechanisms responsible 
for decoherence are not easily individuated and examples 
are provided by two recent experiments which observed 
Rabi oscillations between two circular states of a Rydberg 
atom in a high-Q cavity M , and of two internal states of 
a ^Be''" ion coupled with the vibrations in a trapping 
potential 0. In both cases one observes damped oscilla- 
tions to a steady state in which the population of each of 
the two levels approaches 1/2. A number of candidates 
have been already considered as possible physical sources 
of decoherence in these cases. In the trapped ion case of 
Ref. H , fluctuations of classical parameters such as the 
intensity of the laser beams used to couple internal and 
vibrational degrees of freedom, or the voltage and the fre- 
quency of the trapping potentials have been suggested. 
In the case of the Rydberg atom in a high-Q cavity, dark 
counts of the atomic detectors, dephasing collisions with 
background gas or stray magnetic fields within the cav- 
ity |Q,g| have been proposed as possible sources of deco- 
herence. Despite this, a complete quantitative explana- 
tion of the observed decay rate of the Rabi oscillations 



(see Ref. |7l) in the two experiments is still lacking. The 
only established fact is that, differently from Ref. [Q, in 
both cases, decoherence has a non- dissipative origin. In 
fact, the observed decay of the Rabi oscillations is much 
faster than the energy relaxation rate in these experi- 
mental configurations. Moreover, the fact that in both 
cases the population of each of the two levels asymptot- 
ically approaches 1/2 cannot be explained in terms of 
dissipative mechanisms as the photon leakage out of the 
cavity. 

A different approach to decoherence has been proposed 
in Q , where a model-independent formalism has been de- 
rived to describe decoherence. Here we shall adopt a 
more pragmatic point of view and we shall use this for- 
malism to explain in simple terms both Rabi oscillation 
experiments, even though they are realized in different 
physical situations. The idea underlying the approach 
of Ref. S is the fact that the interaction time, i.e. the 
time interval in which the effective Hamiltonian evolution 
takes place, is a random variable. This randomness can 
have different origins depending on the studied system. 
For example, in the case of the Rydberg atom experi- 
ment M , the interaction time is determined by the tran- 
sit time of the velocity-selected atom through the high-Q 
microwave cavity. This interaction time is random, due 
to fluctuations of the atomic velocities. This random- 
ness implies having random phases e"*^"*/'' in the energy 
eigenstates basis. The experimental results unavoidably 
average over these random phases and this leads to deco- 
herence, i.e., to the decay of off-diagonal matrix elements 
of the density operator in the energy basis. Notice how- 
ever that one would have the same phase fluctuations 
if the Hamiltonian (and therefore the eigenvalues En) 
fluctuates instead of the interaction time. Therefore, as 
we shall see, our approach will give a generalized phase- 
destroying master equation, able to describe many situ- 
ations in which decoherence is associated with random 
phases, originating for example from some frequency or 
interaction time fluctuations. 

Let us consider an initial state p{Qi) and consider the 
case of a random evolution time. The experimentally 
observed state is not described by the usual density ma- 
trix of the whole system p{t), but by its time averaged 
counterpart H 



P{t) 



dt'P{t,t')p{t') 



(1) 



where p{t') — exp{— z_Lt'}p(0) is the usual unitarily 
evolved density operator from the initial state and L . . . — 



[H, . . .]/h. Hence one can write 
p{t) = V{t)p{0) 



(2) 



where V{t) = J^ dt'P{t,t')e-'^*' . In Ref. [|, the func- 
tion P{t, t') has been determined so to satisfy the fol- 
lowing conditions: i) p{t) must be a density operator, 
i.e. it must be self-adjoint, positive-definite, and with 
unit-trace. This leads to the condition that P{t, t') 
must be non-negative and normalized, i.e a probabil- 
ity density in t' so that Eq. (|^) is a completely posi- 
tive mapping, ii) V{t) satisfies the semigroup property 
V{ti +t2) = V{ti)V{t2X, with ti,t2 > 0. These require- 
ments are satisfied by 



V{t) 
P{t,t') 



{l + iLry 



-t/j 



-t'/r 



it'M 



{t/r 



nt/r) 



(3) 
(4) 



so that V{t) and P{t,t') are connected by the so-called 
F- function integral identity [^|lO| . The parameter t char- 
acterizes the strength of the evolution time fluctuations. 
When T ^ 0, P{t, t') -^ 5{t - t') so that p{t) = p{t) and 
V{t) — cxp{— iLt} is the usual unitary evolution. How- 
ever, for finite r, the evolution operator V{t) of Eq. (g) 
describes a decay of the off-diagonal matrix elements in 
the energy representation, whereas the diagonal matrix 
elements remain constant, i.e. the energy is still a con- 
stant of motion. In fact, by differentiating with respect to 
time Eq. (g) and using (||) , one gets the following master 
equation for p{t) 



'p{t) = —\og{l + iLT)p{t) 



(5) 



If one expands the logarithm at second order in r, one 
obtains 



■p{t)^~-[H,-p{t)] 



'[H,[Hrpm 



(6) 



which is the well-known phase-destroying master equa- 
tion [nl| . Hence Eq. (1q) appears as a generalized phase- 
destroying master equation taking into account higher 
order terms in r. Notice, however, that the present ap- 
proach is different from the usual master equation ap- 
proach in the sense that no perturbative and specific sta- 
tistical assumptions are made. 

We now apply this formalism to the two experiments 
of Refs. [||j§|. In the experiment of Ref. [Q, the reso- 
nant interaction between a quantized mode in a high-Q 
microwave cavity (with annihilation operator a) and two 
circular Rydberg states (|e) and \g)) of a Rb atom is 
studied. This interaction is well described by the usual 
Jaynes-Cummings |l3] model, which in the interaction 
picture reads 



H = hnR{\e){g\a+\g){e\a^) 



(7) 



where £7/? is the Rabi frequency. The Rabi oscillations 
describing the exchange of excitations between atom and 



cavity mode are studied by injecting the velocity-selected 
Rydberg atom, prepared in the excited state |e), in the 
high-Q cavity and measuring the population of the lower 
atomic level g, Peg{i) as a function of the interaction 
time i, which is varied by changing the Rydberg atom 
velocity. In the case of vacuum state induced Rabi oscil- 
lations, the decoherence effect is particularly evident and 
the Hamiltonian evolution according to Eq. (0) predicts 



Pea{t) = -{l-cos{2n„t)) 



(8) 



Experimentally instead, damped oscillations are ob- 
served, which are well fitted by 



ptl'd) = J 



^4(1 



-7*, 



1 {2^Bt)) 



(9) 



where the decay time fitting the experimental data is 
7~^ = 40/xsec B and the corresponding Rabi frequency 
is J7fl/27r = 25 Khz. This decay of quantum coherence 
cannot be associated with photon leakage out of the cav- 
ity because the cavity relaxation time is larger (220 /^sec) 
and also because in this case one would have an asymp- 
totic limit P^^P{oo) = 1. The damped behavior of Eq. (|^) 
is instead easily obtained if one applies the approach de- 
scribed above. In fact, from the linearity of Eq. (|^), 
one has that the time averaging procedure is also valid 
for mean values and matrix elements of each subsystem. 
Therefore one has 



Peg{t) 



dt' P{t,t')Peg{t' 



(10) 



Using Eqs. (|2|), (|), (|) and (|), Eq. ^ can be rewritten 
in the same form of Eq. (0) 



where 



Peait) ^- {I -e--^' cos {vt)) 



7 = -log(l + 4f7?jr2) 
V = — arctg (2rifl,T) 

T 



(11) 

(12) 
(13) 



We note that in general the time averaging procedure 
introduces not only a damping of the probability oscilla- 
tions but also a frequency shift. However, if the charac- 
teristic time T is sufficiently small, i.e. VIrt ^ 1, there 
is no phase shift, u ~ 217^^, and 



7 



2f7?,r 



(14) 



The fact that in Ref. M the Rabi oscillation frequency 
essentially coincides with the theoretically expected one, 
suggests that the time r characterizing the fluctuations of 
the interaction time is sufficiently small so that it is rea- 
sonable to use Eq. (H). Using the above values for 7 and 
VLr, one can derive an estimate for r, so to get r ~ 0.5 



/xsec. This estimate is consistent with the assumption 
flRT <C 1 we have made, but, more importantly, it turns 
out to be comparable to the experimental value of the un- 
certainty in the interaction time. In fact, the fluctuations 
of the interaction time are mainly due to the experimen- 
tal uncertainty of the atomic velocity v. In fact, one has 
t — y^w/f , where w is the cavity mode waist. Since 
w = 0.6 cm, the mean velocity is u ~ 300 m/sec and 
the velocity uncertainty is Sv/v = 1% (see Ref. Q), one 
has i = y/nw/v ~ 50 /isec and r ~ (5t = iSv/v — 0.5 
^sec, which is just the estimate we have derived from 
the experimental values. This simple argument supports 
the interpretation that the decoherence observed in [Q is 
essentially due to the randomness of the interaction time. 
Let us now consider the case of the trapped ion ex- 
periment of Ref. p], in which the interaction between 
two internal states (| t) and | |)) of a Be ion and the 
center-of-mass vibrations in the z direction, induced by 
two driving Raman lasers is studied. In the interaction 
picture with respect to the free vibrational and internal 
Hamiltonian, this interaction is described by the follow- 
ing Hamiltonian M 



H ^hn\ t)(i \e^p{i[r]{ae' 
+ H.C. , 



■ a'e * 



)-5t + 



(15) 



where a denotes the annihiliation operator for the vibra- 
tions along the z direction, ojz is the corresponding fre- 
quency and 5 is the detuning between the internal transi- 
tion and the frequency difference between the two Raman 
lasers. The Rabi frequency Jl is proportional to the two 
Raman laser intensities, and 77 is the Lamb-Dicke param- 
eter Wh- When the two Raman lasers are tuned to the 



first blue sideband, i.e. 5 



Hamiltonian (n5[) pre- 



dicts Rabi oscillations between | i, n) and | 1,71 + 1) 
is a vibrational Fock state) with a frequency W 



U U ^l " — -~ u u 






vLliv^) 



(16) 



where L\ is the generalized Laguerre polynomial. These 
Rabi oscillations have been experimentally verified by 
preparing the initial state | i,n), (with n ranging from 
to 16) and measuring the probability Pi{t) as a function 
of the interaction time i, which is varied by changing 
the duration of the Raman laser pulses. Again, as in 
the cavity QED experiment of y] , the experimental Rabi 
oscillations are damped and well fitted by ||,R] 



nKi)-2(i 



-fr.t , 



's(2aii)) 



(17) 



where the measured oscillation frequencies r2„ are in 
very good agreement with the theoretical prediction ( |lq ) 
corresponding to the measured Lamb-Dicke parameter 
•q — 0.202 Q. As concerns the decay rates 7„, the exper- 
imental values are fitted in g by 



ln=lo{n + lf 



(18) 



where 70 = 11.9 Khz. This power-law scaling has been 
investigated in Refs. ||l^,Q, but a clear explanation of 
this behavior of the decay rates is still lacking. On the 
contrary, the scaling law (|l^) can be accounted for in 
the previous formalism if we consider the small r limit of 
Eq. (Q), which is again suggested by the fact that the 
experimental and theoretical predictions for the frequen- 
cies Vtn agree. In fact, the n-dependence of the theoretical 
prediction of Eq. (n6h for 77 = 0.202 is well approximated, 
within 10 %, by the power law dependence 



n^~Vto{n + l) 



0.35 



(19) 



so that, using Eq. ( |l4| ) with Q.^ replaced by rj„, one 
has immediately the power law dependence (71 -I- l)*^'^ of 
Eq. (|lj) . The value of the parameter r can be obtained 
by matching the values corresponding to 7i = 0, and using 
Eq. (|l|), that is r = 7o/20g ~ 1.5 • 10"^ sec, where we 
have used the experimental value rio/27r = 94 Khz. How- 
ever, this value of the parameter r cannot be explained 
in terms of some interaction time uncertainty, such as 
the time jitter of the Raman laser pulses, which is ex- 
perimentally found to be much smaller |15| . In this case, 
instead, the observed decoherence can be attributed, as 
already suggested in |0,|l3|,0 , to the fluctuation of the 
Raman laser intensities, yielding a fluctuating Rabi fre- 
quency parameter U,(t) of the Hamiltonian (15). In this 
case the evolution is driven by a fluctuating Hamiltonian 
H{t) = nn{t)H, where H = H/Q in Eq. (P), so that 



p{t) = exp -iL / d^niO p(0) 



_ -2LA{t) 



p(0) (20) 



where L = [H, . . .]/?i and we have defined the positive 
dimensionless random variable A(t) = L d^^{^), which 
is proportional to the pulse area. It is now easy to un- 
derstand that the physical situation is analogous to that 
characterized by a random interaction time considered 
above, with L replaced by L and t' by A{t). One has 
again phase fiuctuations in the energy basis representa- 
tion and, in analogy with Eq. (|l|), one considers an aver- 
aged density matrix 



m 



dAP{t, A)e-'-^'^p(O) 



(21) 



Imposing again that p{t) must be a density operator and 
the semigroup property, one finds results analogous to 
Eqs. d) and g) 



-t/r 



V{t) = f 1 + iLnrj 



Pit, A) 



nr 



r(t/7 



(22) 
(23) 



where, the parameter Q of Eq. dla) plays now the role 
of a mean Rabi frequency. In fact, consistently with the 
probability distribution of Eq. (p3|), one has il = {A)/t. 



The scaling time r characterizes in this case the strength 
of the pulse area fluctuations, since from Eq. (|2J^), one 
has 0-2 (A) = {A^) - (A) 2 = 9?-tT. The estimated value 
of T is reasonable since it corresponds to a fractional er- 
ror of the pulse area yJa'^{A)/{A) — \/t jt of 10% for a 
pulse duration oil = \ fisec, and which is decreasing for 
increasing pulse durations. 

The present analysis shows many similarities with that 
of Ref. n^ which also tries to explain the decay of the 
Rabi oscillations in the ion trap experiments of |q] in 
terms of laser intensity fluctuations. The authors of 
Ref. 1 13 1 in fact use a phase destroying master equa- 
tion coinciding with the second-order expansion (m of 
our generalized master equation of Eq. (||) (see Eq. (16) 
of Ref. (l3| with the identifications G ^> H/h and T ^^ t) 
and moreover derive the same numerical estimate for the 
pulse area fluctuation strength F ^-> r. Despite this simi- 
larities, they do not recover the scaling (nsh of the decay 
rates 7„ only because they do not use the general ex- 
pression of the Rabi frequency (M), (and which is well 
approximated by the power law (R%) but its Lamb-Dicke 
limit fin = flo{n+ 1)°'^, which is valid only when 77 <IC 1. 
There is however another, more fundamental, difference 
between our approach and that of Ref. |Q. They as- 
sume from the beginning that the laser intensity fluc- 
tuations have a white and gaussian character, while we 
make no a priori assumption on the statistical proper- 
ties of the pulse area A. We derive these properties, i.e. 
the probability distribution (E3) , only from the very gen- 
eral semigroup condition, and it is interesting to note 
that this condition yields a gaussian probability distri- 
bution for the pulse area only as a limiting case. In fact, 
from Eq. ( p^ ) one can see that P(t, A) tends to become 
a gaussian with the same mean value Qt and the same 
variance fl^rt only in the large time limit t/r ^ 1. This 
non-gaussian character of P{t,A) can be traced back to 
the fact that P(t, A) must be definite and normalized in 
the interval < A < -f-oo and not in —00 < A < -|-cx). 
Notice that a,t t — t, Eq. 



form P{t = T,A) 



'Al 



23) assumes the exponential 
/rJr. Only at large times 



t the random variable A becomes the sum of many in- 
dependent contributions and assumes the gaussian form. 
Due to the non-gaussian nature of the random variable 
A^ we find that the more generally valid phase-destroying 
master equation is given by Eq. (m (with L replaced by 
f2L), whose predictions significantly depart from its sec- 
ond order expansion (g) (corresponding to the gaussian 
limit) as soon as r becomes comparable with the typi- 
cal timescale of the system under study, which, in the 
present case, is the inverse of the Rabi frequency. 

In conclusion, we have presented a model-independent 
theory for non-dissipative decoherence, able to provide a 
simple and unified description of the same decoherence 
phenomenon observed in two Rabi oscillations experi- 
ments which were performed under different situations. 
A simple way to test experimentally our prediction is to 
check that the dependence of the decay rate as a function 



of the Rabi frequency is given by Eq. ( p^ ) . One should 
observe a transition from a quadratic dependence to a 
logarithmic dependence, increasing the value of the Rabi 
frequency, or of t. 

Discussions with J.M. Raimond and D. Wineland are 
greatly acknowledged. This work has been supported by 
MURST under the "Cofinanziamento" . 



[1] W.H. Zurek, Pliys. Today 44(10), 36 (1991), and refer- 
ences therein. 

[2] M. Brune et al, Phys. Rev. Lett. 77, 4887 (1996). 

[3] D. Vitali, P. Tombesi, G.J. Milburn, Phys. Rev. Lett. 79 
2442 (1997); Phys. Rev. A 57 4930 (1998). 

[4] M. Brune et al, Phys. Rev. Lett. 76, 1800 (1996). 

[5] D.M. Meekhof et al, Phys. Rev. Lett. 76, 1796 (1996). 

[6] J.M. Raimond, private communication. 

[7] D.J. Wineland et al, J. Res. Natl. Inst. Stand. Technol. 
103, 259 (1998). 

[8] R. Bonifacio, N uovo Cimento 114 B, 473 (1999); LANL 
e-print archive quant-ph/990106S ; in Mysteries, Puz- 
zles and Paradoxes in Quantum Mechanics, edited by R. 
Bonifacio, AIP, Woodbury, 1999, pag. 122. 

[9] I.S. Gradshteyn and I.M. Ryzhik, Table of Integrals and 
Series, Academic, Orlando, 1980, pag. 317. 
[10] In Ref. S a more general expression for P{t, t') and V{t), 
depending on two parameters ri and T2 is derived. We 
choose T\ = T2 = T because in the experiments considered 
here, the effective interaction time (t') — iri/r2 |8| has to 
coincide with the "laboratory time" t, implying therefore 
n = T2. 
[11] D.F. Walls and G.J. Milburn, Quantum Optics, Springer, 

Berlin, 1996. 
[12] E.T. Jaynes and F.W. Cummings, Proc. IEEE 51, 89 

(1963). 
[13] S. Schneider and G.J. Milburn, Phys. Rev. A 57, 3748 

(1998). 
[14] M. Murao, P.L. Knight, Phys. Rev. A 58, 663 (1998). 
[15] D. Wineland, private communication. 



